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Abstract 

We analyze the time-dependent solutions of the pseudo-differential Levy- 
Schrodinger wave equation in the free case, and we compare them with the 
associated Levy processes. We list the principal laws used to describe the 
time evolutions of both the Levy process densities, and the Levy-Schrodinger 
wave packets. To have self-adjoint generators and unitary evolutions we will 
consider only absolutely continuous, infinitely divisible Levy noises with laws 
symmetric under change of sign of the independent variable. We then show 
several examples of the characteristic behavior of the Levy-Schrodinger wave 
packets, and in particular of the bi-modality arising in their evolutions: a 
feature at variance with the typical diffusive uni-modality of both the Levy 
process densities, and the usual Schrodinger wave functions. 

PACS: 05.40.Fb, 02.50.Ga, 02.50.Ey 
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1 Introduction and notations 

In a recent paper pQ it has been shown how to extend the well known relation 
between the Wiener process and the Schrodinger equation [2],[3jlH[5] to other suitable 
Levy process. This idea - discussed elsewhere only in the stable case P, [7] - leads 
to a L-S (Levy-Schrodinger) equation containing additional integral terms which 
take into account the possible jumping part of the background noise. In fact, the 
infinitesimal generator of the Brownian semigroup (the Laplacian) being substituted 
by the more general generator of a Levy semigroup, we get an integro-differential 
operator with both a continuous (differential and Gaussian) and a jumping (integral, 
non Gaussian) part. These ideas have already been discussed in the framework of 
stochastic mechanics [21 [5] and are considered as a model for systems more general 
than just the usual quantum mechanics: a true dynamical theory of Levy processes 
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that can be applied to several physical problems jH]- The aim of this paper is now to 
show a number of explicit examples of wave packets solutions of L-S free equations. 

In recent years we have witnessed a considerable growth of interest in non Gaus- 
sian stochastic processes - and in particular into Levy processes - from statistical 
mechanics to mathematical finance. In the physical field, however, the research 
scope is presently rather confined to the stable processes and to the corresponding 
fractional calculus [SI El E], while in the financial domain a vastly more general type 
of processes is at present in use. Here we suggest that a Levy stochastic mechanics 
should be considered as a dynamical theory of the entire gamut of the infinitely divis- 
ible processes with time reversal invariance, and that the horizon of its applications 
should be widened even to cases different from the quantum systems. 

This approach has several advantages: first of all the use of general infinitely 
divisible processes lends the possibility of having realistic, finite variances. Second, 
the presence of a Gaussian component and the wide spectrum of decay velocities 
of the increment densities will give the possibility of having models with differences 
from the usual Brownian (and usual quantum mechanical, Schrodinger) small 
as we want. Last but not least, there are examples of non stable Levy processes 
which are connected with the simplest form of the quantum, relativistic Schrodinger 
equation: an important link that was missing in the original Nelson model. This final 
remark, on the other hand, shows that the present inquiry is not only justified by 
the a desire of formal generalization, but is required by the need to attain physically 
meaningful cases that otherwise would not be contemplated in the narrower precinct 
of the stable laws. 

In this paper we will show practical examples for the behavior of the evolving 
wave packet solutions of particular kinds of (non Wiener) L-S equations, and we 
will put in evidence their characteristics. In particular the bi-modality arising in 
many of these these evolutions which has a correspondence neither in the the process 
diffusions, nor in the usual Schrodinger wave functions: an effect which has already 
been observed only in confined Levy flights [10J. This is coherent with the usual 
stochastic mechanics scheme, in so far as in this theory the Schrodinger equation 
is recovered by introducing a kind of interaction modeled by means of a quantum 
potential [2J [5]. In the following exposition laws and processes will always be one 
dimensional. An extensive analysis of the topics discussed in this first chapter is 
available in the two monographs [UJ and [12] , while a short introduction can be 
found in [13] . 

In the present paper the law of a rv (random variable) X is characterized either 
by its pdf (probability density function) /, when - as it is generally supposed - the 
law is ac (absolutely continuous), or by its chf (characteristic function) (p with the 
usual reciprocity relations 



When the laws are not ac we sometimes will use the Dirac delta notation: the 
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symbol S xo (x) = S(x — Xq) will then represent a law degenerated in x and only 
formally it will act as a pdf. The symbol S(x) will also be used instead of 5q(x). 
In order to have background noises with generators self-adjoint in L 2 - an essential 
requirement for our purposes - we will consider only symmetric laws, namely we 
will require 

f(~ x ) = f( x ) > <p(-v) = Viu) 

so that the chf <£> will also be real. This also means that, when it exists, the 
expectation vanishes E [X] = 0, namely the law is also centered. See the Appendix [Al 
for further details about our notations. 

Since we will restrict our analysis to background noises driven by Levy processes, 
we will be interested almost exclusively in id (infinitely divisibld3, for details see [TTj 
[T2j IT5] ) laws with a Levy triplet C = (a,/3,v). Here our Levy measures v will 
always be supposed to have a density: v{dy) = £(y) dy; when this does not happen 
we will often use the Dirac delta notation. As a consequence the Levy triplet will 
be rather specified as C = (a,(3,£). The Ich (logarithmic characteristic) of our id 
laws rj = hi if, with ip = e v , will then satisfy the Levy-Khintchin formula 

V (u) =iau-\ (3 2 u 2 + / [e my - 1 - tuy I D (y)] £(y) dy (2) 
1 Jy^o 

where D = {y : \y\ < 1}. The prescription of the integral around the origin is 
essential only when - as usually may happen - the Levy measure shows a singularity 
in y = 0. When the law is dimensionless (see Appendix |A|) then also a, (3, I and y 
are so; on the other hand, if the law has the dimensions of a length, then a, /3, y are 
lengths, while i is the reciprocal of a length. In particular when the law is symmetric 
we have 

or = 0, £(-x)=£(x) 
so that the Levy-Khintchin formula will be reduced to the symmetric real expression 

V (u) = ~/3 2 u 2 + [ (cosuy-l)£(y)dy (3) 

and hence the chf p will not only be real, but also non negative: f(u) > 0. 

The Markov processes dealt with in this paper are stationary, independent in- 
crements processes and are then defined by means of the chf ip M ' r of their At- 
increments, where r is a dimensional, time scale parameter. Here too we can intro- 
duce a dimensionless formulation through the coordinate t/r, and to simplify the 
notation we can continue to use the same symbol t for this dimensionless time. In 
this case the stationary chf will be reduced to <p At , and the dimensional formulation 

X A law (p is said to be id if for every n it exists a chf tp n such that tp = <^™; on the other hand 
tp is said to be stable when for every c > it is always possible to find a > and b G R such that 
e lbu ip(au) = [<p(u)] c . Every stable law is also id. See also Appendix [S] for further details about 
stable laws. 
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will be recovered by simple substitution of t/r to t. A stochastically continuous 
process with stationary and independent increments is called a Levy process when 
X(0) = 0, P-a.s., but this paper will mostly be about the same kind of processes for 
arbitrary initial conditions X(0) = X , P-a.s. with law fo(x) and fo(u) = e v °( u \ All 
these processes, independently from their initial conditions, will share both the same 
differential equations (whether SDE's, or PDE's) and the same transition pdf's 

f X (t)(x\X(s) = y) = p(x,t\y,s). 

To avoid confusion we will then adopt different notations for their respective marginal 
pdf's: for a Levy process (namely with Xq = initial condition) we will write 

fx(t) (x) =q(x,t), ip X (t) (u) =x(u,t) 

with q(x, 0) = 5(x) and x{ x i 0) = 1> while for the general stationary and independent 
increments process (with arbitrary initial condition Xq) we will write 

fx(t) (x) = p(x, t) , ifx{t) (u) = <t>(u,t) 

with p(x, 0) = fo(x) and <f>(x, 0) = (po(x). It is then easy to show that 

P(x,t\y, s) = q(x-y,t-s). (4) 

The infinitesimal generator A = i](d) (here d stands for the derivation with 
respect to the variable of a test function v) of the semigroup of a Levy process will 
be a pseudo-differential operator with symbol r\ [U [12] , namely from (j2J) 

[Av](x) = [r)(d)v](x) = —j= \ e iux 7](u)v(u) du 

V27T 7_no 



+ / [v(x + y)- v(x) - yI D (y) d x v(x)} £(y) dy 

where v denotes the FT (Fourier transform) of the test function v with the reci- 
procity relations: 

v(u) = j v(x)e~ lux dx, v(x) = —= / v(u)e lux du 

V27T J-oo V27T J-oo 



The generator A will be self-adjoint in L 2 (R,dx) when the law is symmetric, and 
in this case it reduces to 

[Av){x] = ^ d 2 x v{x) + / [v{x + y)- v{x)} £{y) dy (6) 
1 Jy^o 
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Table 1: List of the essential properties of a few basic, dimensionless laws discussed 
in this paper: degenerate (Dirac) D, normal (Gauss) 91, Cauchy £, Laplace £, 
uniform it, and doubly degenerate in +1,-1 (symmetric Bernoulli) D%. 



so that it is determined by the two essential elements of our Levy triplet, namely f3 
and £. Given the process stationarity, in a dimensionless formulation the transition 
law degenerate in x = at t = will have as chf x — V 9 * = etv an d as pdf 

i p+oo i r+oo 

q(x, t) = — X (u, t)e~ iu * du=— tp(u)*e*" du. (7) 

J_ 00 27T J_ OQ 

This transition law plays an important role in the evolution of an arbitrary initial 
law /o, <£>o : the process chf will indeed be now <p(u,t) = x(% t)fo{ u )-i an d the 
corresponding pdf will be calculated from 

i r>+oo 

p(x,t) = [q(t) * f ](x) = — / <P(u,t)e-^du 

In 



— oo 



namely either as a convolution of the transition and the initial pdf's, or by inverting 
the chf (f> of the process. This pdf will also be a solution of the evolution pseudo- 
differential equation [TJ [T2] 

d t p = v(d)p, p(x, 0) = f (x) (8) 

which from (jSJ) takes the integro-differential form 

(3 2 

d t p(x, t) = a d x p(x, t) + — d%p(x, t) 

+ I \p( x + y> *) _ P( x > *) ~ y J D(y) d x p(x, t)} £(y) dy 

Jy^O 
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Table 2: List of the generators of the Levy processes associated to some of the non 
degenerate, id, dimensionless laws of Table Q] 

and for a centered, symmetric noise from reduces to 

d t p{x,t) = ^-d 2 x p{x,t) + [p{x + y,t)-p{x,t)}e{y)dy. (9) 

We finally remember that, since (EJ) and ([9]) are given in terms of process pdf 8, 
this equations are supposed to hold only for ac processes. We are then required to 
point out which Levy processes have densities. To answer - at least partially - this 
question we then recall that [11] any non-degenerate, sd (self-decomposable^], for 
details see [TTJ [121 Q2]) distribution is ac. On the other hand such a property also 
extends to the corresponding processes for every t. In fact [11] if X(t) is a sd process 
also its pdf at every t is sd, and hence X(t) is ac for every t. As a consequence 
we can always explicitly write down the evolution equations ([§]) in terms of the 
process pdf's at least for the sd case. We remark, however, that there are also non 
sd processes which are ac: the ac compound Poisson processes of Appendix |B1 are 
an example in point. 

We listed in the Tabled] the properties of a few basic, symmetric, dimensionless 
laws: degenerate (Dirac) D, normal (Gauss) 91, Cauchy £, Laplace £, uniform it, 
and doubly degenerate in +1,-1 (symmetric Bernoulli) D\. The uniform law pdf 
is given by means of the Heaviside functions Q(x). These laws are also relevant 
particular cases of the families that we will introduce in the Section [2] Remark that 
in the Table [1] there is no value for the expectation of £ because it does not exist 
(£ is centered on the median), and no values for the Levy triplet of il and T>i since 
these are not id laws. Moreover in general our laws are not necessarily standard. 
The form of the simplest generators corresponding to our Levy processes is finally 
shown in the Table [2] 

The paper is organized as follows: in the Chapter [2] we recall the essential prop- 
erties of the law families of our interest; then in the Chapter [3] the L-S equation 
is introduced with its connections to the Levy processes. Finally in Chapter H] our 

2 A law ip(u) is sd when for every a £ (0, 1) we can always find another chf tp a (u) such that 
<p(u) = (p(au)(p a (u). Every stable law is also sd; every sd law is also id. 
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Figure 1: Graphical synthesis of the relations among the families of laws discussed 
in the Section |2j The uniform it is our unique example beyond the pale of the id 
laws, while the laws of the (simple) Poisson family 23(A) are id but not sd. Notable 
cases (9T, £, (£) within the sd families are put in evidence; the Cauchy € law lies at 
the intersection of the stable 6(A) and Student T(A) families. 

examples are elaborated and in Chapter [5] the results are collected and discussed. A 
few technical details are collected in the Appendices in order to avoid to excessively 
burden the text. 

2 Families of id laws 

We will introduce here the principal families of id laws considered in this paper. For 
a graphical synthesis of the relations among them see Figure [TJ Please remark that 
this synthesis is particularly simple because we limit ourselves here to dimensionless 
laws (see Appendix |AJ: this produces one-parameter families that can be easily 
represented in our scheme. In the Table [3] are then listed the properties of the 
principal families of dimensionless, sd laws that will be discussed. The " . . . " symbol 
in this table means either that we do not have an elementary formulation for the 
entry, or that there are no particular values of A to be put in evidence. K u , B and 
T respectively are the modified Bessel functions of the second kind, and the Euler 
Beta and Gamma functions, while H\ stands for the Fox if -functions representing 
the pdf of stable laws [TT] • From Table [1] and Table [3J we can on the other hand 
immediately see that 6(1) = = €, 6(2) = 91 and 23(25(1) = £, as also put in 
evidence in the Figure [TJ The behaviors of a few /c/i's of id laws are finally displayed 
and compared in the Figure [2j it could be seen there that all the Ich's of the sd 
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Figure 2: The Ich —rj(u) of some basic dimensionless laws from Tabled plus that 
of a compound Poisson ^(A,^) with normal component laws (see Appendix IB]) . 
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Table 3: Properties of our principal families of sd, dimensionless laws: the sta- 
ble ©(A), the Variance-Gamma QJ(5(A), the Student T(A) and the relativistic qm 
(quantum mechanics) 9t(A). 
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laws considered in this paper diverge at infinity with velocities ranging from u 2 to 
logw, while the unique not diverging Ich characterizes one of our non sd examples: 
the compound Poisson ^3(A, OT) with normal component laws. This also gives an 
intuitive idea of how much the behavior of a law - in so far as we are concerned, 
for instance, with its jumping properties - differs from that of the P-a.s. continuous 
Gaussian case. 



2.1 The stable laws 6(A) 

This is the more widely studied family of id laws, albeit among them only the normal 
(3(2) = enjoys a finite variance. But for the 9T, the £ and precious few other cases 
the pdfs of the stable laws exist only in the form of Fox if-functions [Hj. To see in 
what sense these laws are stable we must for a moment reintroduce the dimensional 
parameter a: we then have a larger family & a (X) with two parameters, < A < 2 
and a > 0, and 

<p( u ) = e - flA H A /\ ( 10 ) 

Now, for a given fixed A, the family <5 a (A) with a > is closed under convolution, 
as can be easily seen from (fTUI) . For instance the families of the normal 9t a = 91(a 2 ) 
and Cauchy £ a laws are closed under convolution since DT(a^)* ^(a 2 ,) = 9T(a 2 + a 2 ,) 
and (£„, * (T a2 — (£ ai _|_ a2 . Stability however means more: the families (3 a (A) for a 
given A are types of laws, in the sense that a law of the family differs from another 
just by a re-scaling (centering is not necessary here because our laws already are 
centered; for details see Appendix [AJ , the parameter a being indeed nothing else 
than a space scale parameter. This has far reaching consequences. In particular it 
is at the root of the well known fact that the stable Levy processes are self -similar: 
a property not extended to other, non stable Levy processes |13j. The generators of 
the stable Levy processes are 



— 1 f v(x + y) — v(x) 

2Ar(-A) cosf Jy^o \y\ 



woo = „^, " ui* *y 0<A<2 > A ^ 



while for A = 1 (<£ law) and A = 2 (9? law) they are listed in the Table [2 



2.2 The Variance— Gamma laws 070(A) 

The Variance-Gamma laws owe their name to the fact that they can be seen as 
normal variance-mean mixture^ where the mixing density is a gamma distribution. 

3 A normal variance- mean mixture, with mixing probability density g, is the law of a random 
variable Y of the form Y = a + jSV + a\/VX where a and /3 are real numbers and a > 0. 
The random variables X and V are independent; A is a normal standard, and V has a pdf g 
with support on the positive half-axis. The conditional distribution of Y given V is then a normal 
distribution with mean a + /3V and variance a 2 V. A normal variance-mean mixture can be thought 
of as the distribution of a certain quantity in an inhomogeneous population consisting of many 
different normally distributed sub-populations. 
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It is apparent moreover from the Table [3] that 030(A) is closed under convolution in 
the sense that 23(5(Ai)*93<5(A2) = 23(5(Ai+A 2 ). That notwithstanding, however, the 
Variance-Gamma laws are not stable. To see that let us reintroduce the dimensional 
scale parameter a to have the enlarged family Q3(5 a (A): 

/ 1 \ A 

(p(u) 



1 + a 2 u 2 / 

Now every sub-family with a given, fixed a is closed under convolution, but at 
variance with the stable case the parameter describing the sub-family is A, rather 
than a. As a consequence the closed subfamilies do not constitute types of laws 
differing only by a rescaling, and hence the laws are not stable. The pdf's of the 
Variance-Gamma laws can be given in particular instances as finite combinations 
of elementary functions. By generalizing the quoted example of the Laplace law 
23(5(1) = £, when A = n + 1 with n = 0, 1, ... we have for the dimensionless pdfs 

fix) = f( 2n - k ) ( 2 N)^ N = J—Ld (Id) 

k=0 v 7 

where 6 n (x) are reverse Bessel polynomials [15]. All our dimensionless 23(5 (A) laws 
are endowed with expectations (which vanish by symmetry) and finite variances 2A. 
The generator of the corresponding Levy process is 

[Av](z) = \f + e -\vl d y A>0 



ly^o \y\ 

which coincides with that of £ (see Table [2]) for A = 1. 
2.3 The Student laws 1(A) 

But for the Cauchy € case, the laws of the Student family (even enlarged by means of 
the scale parameter a) are not stable, and T(A) itself is not closed under convolution: 
convolutions of Student laws are not Student laws. As can be seen from Table [3] 
the Variance-Gamma and the Student families enjoy a sort of duality since their 
pdf's and c/i/'s are essentially exchanged. This has been discussed at length in a 
few recent papers [T7J HH [19]. Remark that to put in evidence this correspondence 
we have chosen the Student laws of T(A) without introducing the usual parametric 
scaling x 2 / A of its variable that would have put equal to A/(A — 2) all their variances 
for A > 2. In particular this means that for A — > +oo we will not get a standard 
9T law, as also shown in the Figure [TJ The following remarks are however virtually 
untouched by this choice. While the pdfs and c/i/'s of the Student laws are known, 
differently from the Variance-Gamma laws their Levy measures and generators have 
not a known general expression. However we can give them in particular instances. 
For example when A = 2n + 1 with n = 0, 1, . . . the chf becomes 

.. A n \(2n-h)\ tn . .. k _, u , n!2 m e-H n .. 

^ = E (2n) i (TI - k )\k\ (2N) e ll = i^r W) ' 
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where 9 n are again reverse Bessel polynomials [15] . Of course T(l) = £ is the well 
known Cauchy (stable) case, while for T(3) we have 

f( x ) = l(l^) > ^«) = (1 + H)e- N . 

and it can be shown [17] in this case that the Levy measure is 

1 — I a; I (sin \x\ ci \x\ — cos \x\ si \x\) 

e{x) = n 

7TX Z 

where the sine and the cosine integral functions are 

[ + °°smy . [ + °°cosy 
six = — / dy , ci x — — / dy 

J x y Jx y 

The existence of the moments of the X(A) laws depends on the value of the parameter 
A: the n th moment exists if n < A. In particular the expectation exists (and vanishes) 
for A > 1, while the variance exists finite for A > 2 and its value is (A — 2) _1 . The 
generator of the Levy process can finally be explicitly given for T(3) from (ITTj) 

MM = a f k* + „) - , W ] i-\v\(™\y\«\»\ -<™\»\s'\y\) dy . 

: <y 

2.4 The compound Poisson laws 9T CT *^3(A,i}) 

The compound Poisson laws Sjo * *}3(A,ij) are not sd, but they are nevertheless id; 
they are also ac when fj is ac (for details and notations see Appendix [B]). In the 
following examples we will take into account the dimensional parameters a of the 
component laws. Consider now the case Sj = OTo-: then £o(x) = and (3o = a so 
that the Levy triplet of 91^ * ^}3(A, Sy) is C = (0 , a , A h) and the generator is 

<j 2 r 

[Au](:r) = — d%.v(x) + A / [v(x + y) — v(x)]h(y) dy . 
2 // •■ o 

When in particular also fj = 9t a , then the Levy triplet of 9t CT * *P(A, 9t a ) is 



£ = I 0, a, A 6 



-a; 2 /2a 2 



V27ra 2 

and we get a law with the following pdf and Ich 



OC 



^ fe! y/2ir{ka 2 + a 2 ) 2 
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namely a Poisson mixture of centered normal laws 91(A;a 2 + a 2 ). The self-adjoint 
generator then is 



0-2 r+oc e -V 2 /2a 2 

[Av\(x) = — d 2 x v{x) + \ J [v(x + y) - v(x))-j== dy . 

and we could look at it as to a Poisson correction to the Wiener generator, the rel- 
ative weight of these two independent components being ruled by the ratio between 
A and a 2 . 

As another example of ac compound Poisson law let us suppose instead that 
■fto = 9^<t again, but that = D a (see Appendix |B|) . so that the Levy triplet of 
m a * <p(A, © ) now is 

£ =(0 ,a,A^ /a) + Ma;/a) 



2a 



while its pdf and Ich are 



~ \k 1 JL /fc\ e -[x-(fc-2j)a] 2 /2a 2 



fc! 2 fc ^ Vj7 \Z27ra2 



7/(w) = A(cos aw — 1) 



Here the law is again a mixture of normal laws 91(na,<7 2 ), n = 0,±1,... which 
however are now centered around integer multiples of a. The generator finally is 

a 2 2 u(x + a) - 2u(x) + - a) 
[Av](a;) = y %v(x) + A^ ^ ^ ^ '- 

because the integral jump term reduces itself to a finite difference term. 



2.5 The relativistic qm laws £H(A) 

The family of the relativistic qm (quantum mechanics) laws on the other hand is a 
particular case of the well known (centered and symmetric) Generalized-Hyperbolic 
family (see for example [17] and references quoted therein): in fact we have 9^(A) = 
&$) (— |, 1, A), as can be seen by direct inspection of their pdfs and chfs. Remark 
as a consequence that these are not simple Hyperbolic laws that constitute the 
different particular sub-family <&Sj(l,l, A). The name follows from the fact that 
- for a suitable identification of the parameters A and a by means of the particle 
mass m, the velocity of light c and the Planck constant h - its pseudo-differential 
generator 

A = rj(d x ) = mc 2 — ^m 2 c A — c 2 h 2 d 2 

coincides with the Hamiltonian operator of the simplest form of a free relativistic 
Schrodinger equation [TJ [12], and hence its corresponding L-S equation exactly co- 
incides with this free relativistic Schrodinger equation (see Appendix ID.5j) . £K(A) is 
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closed under convolution, as can be seen from the form of the chfs, but the laws 
are not stable for the same reasons as the Variance-Gamma: the parameter A is 
not a scale parameter. The pdf's and c/i/'s are explicitly known (see Table [3]), and 
all their moments exist: the odd moments (in particular the expectation) vanish by 
symmetry, while the even moments are always finite and its variance is A. Since the 
Levy measure is explicitly known (see Table [3D the Levy dimensionless generator 
also takes the form 



where K\ is a modified Bessel function. 

3 The Levy— Schrodinger equation 

To keep the notations as simple as possible also in this chapter the laws and the 
time coordinate will again be supposed dimensionless. It has been shown in [T] that 
the evolution equation (Q of a centered, symmetric Levy process can be formally 
turned into a L-S equation: in fact the pseudo-differential generator t]{9) of our 
processes is a self-adjoint operator in L 2 and hence can correctly play the role of 
a hamiltonian. We summarize in the following the formal steps leading to the L-S 
equation (for further details see [1]); this will also establish the notation for the 
subsequent sections. 

Take as background noise a centered, symmetric, id law with /, ip = e v , £ = 
(0, j3, £) and a symmetric £ so that (|3J) holds 



remember that since r\ is real and symmetric, ip too will be real, symmetric and non 
negative (ip > 0). Define then the transition chf xi u it) = V 9< ( M ) an d the reduced 
transition pdf 



of the corresponding Levy process, and take an initial law / , (po = e Vo : the chf and 
the pdf of the process will be 




dy. 





<p(u,t) 



X{u,t)<p (u) 




(12) 
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There are hence two ways to calculate p(x,t): either as p = q * fo, or by inverting 
the chf (ft = x^Po- As a matter of fact these two ways give the same result, but - 
depending on the specific problem - one can be easier to calculate than the other. 
The pdf p(x,t) of the previous step must also be a solution of the (dimensionless) 
evolution equation 

/3 2 f 

d t p(x,t) = —d 2 x p(x,t) + \p{x + y,t)-p(x,t)]£{y)dy (13) 

and in principle we could find p also by directly solving this equation. 

We pass then to the L-S propagators by means of the formal substitution t — > it: 

j(u,t) = X (u,it) = <p u (u) = e lt ^ u \ g(x,t) = q(x,it) 

so that g and 7 will still verify the same reciprocity relations (J7J) of q and x 

g(x, t) = q(x, it) = — / x(% H)e~ iux du = — 7 («, t)e~ iux du. 



—00 j— 00 



Remark that if the law of the background noise is centered, symmetric and id then 77 
is real, symmetric and positive and hence we always have I7I = 1. This implies first 
that 7 is not normalizable in L 2 , and hence that also g is not normalizable in L 2 . 
This is not surprising since, as it is well known, the propagators are not supposed to 
be normalizable io/'s. On the other hand, as we will see later, this also entails that 
an initial normalized wf will stay normalized all along its evolution. We choose now 
an initial L-S wf: to compare the evolutions of the wfs with that of the process 
pdf's, we will start - whwnever we can - with a law / , (po = e Vo and with a wf ipo 
such that 1 2 = fo, namely 

where So is an arbitrary, dimensionless, real function. In this way we are also sure 
that ipo £ -^ 2 (B0; and that HV'oll 2 = 1- As a matter of fact we could also characterize 
our initial state through the wf FT ipo{u) which exists because ipo £ L 2 (M). It is 
possible to show that (po and ipo must satisfy the following relation 

^0 = i>o * io 

which is simply the dual of \ipo\ 2 — fo- The initial wf can be simplified by choosing 
fo and (po centered and symmetric, with So = 0. In this way we will have real tpo 
and ipo, with 

Mx) = \ r W)> (14) 

so that the following relation will always be satisfied 

¥0 = ^0*^0- (15) 



N CUFARO PETRONI: Levy-Schrodinger wave packets 



15 



Now the L-S wfs will obey the following evolution scheme 



4>(u,t) = -f(u,t)ip {u) 



^OM) = [g(t) * ip ](x) = / g(x - y,t)i/j (y)dy 



J — oo 



ij)(x,t) = -= / ij(u,t)e iux du 




(16) 



Here we can see the relevance of having |7| 2 = 1 (namely of having a centered, 
symmetric background Levy noise, and hence a self-adjoint generator): we have 
indeed that \^(t)\ 2 = |7| 2 |^o| 2 = |V>o| 2 , so that if ||^ || 2 = 1 then also ||^(t)|| 2 = 1, 
and as a consequence (by Parseval and Plancherel theorems) ||^>(£)|| 2 = 1 at every 
t. In other words we can say that the non normalizability of the propagator is the 
counterpart of the unitarity of the L-S evolution. Finally the w/'s ip(x, t) introduced 
in the previous steps must satisfy the free L-S equation 



4 Processes and wave packets 

We will give now several examples of L-S wfs compared with the corresponding 
purely Levy evolutions. We classify these examples first by choosing the laws of the 
background noises: this will be done by picking up the id laws that allow a reasonable 
knowledge of both the transition pdf of the Levy process, and the L-S propagator. 
Besides the usual Wiener case (that will be considered just to show the way) this will 
indeed allow us to calculate the evolutions by means of integrations, without being 
obliged to solve pseudo-differential equations. The equation will be used instead - 
when it is possible - as a check on the solutions found from transition pdfs and 
propagators. We will compare then the typical evolutions of the Levy process pdfs, 
and of the wfs solutions of a free L-S equation: for details, notations and formulas 
about both the initial laws and wfs, and the transition pdfs and propagators we 
will make due references to Appendix O and to Appendix [Dl Remark also that in 
the following we will reintroduce the dimensional parameters a, b and r. 

4.1 Gauss 

Take a Wiener process with transition law ( HHl) : for a normal initial law ( |35i) 9Tb we 
have 



so that the evolution is always Gaussian Dl(2Dt + 6 2 ): it starts with a non degenerate 
normal distribution of variance b 2 and then widens as the usual diffusions do with 
variance 2Dt + b 2 . The L-S evolution of the w/'s on the other hand is here the usual 




(17) 



</>(M) = X(u,t)ip (u) = e 



(2Dt+b 2 )u 2 /2) 
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quantum mechanical one: take first as initial wf the Gaussian (1361) : we then have 
as wave packets 

= 7(M)i(«)=}fe-^« ffl )" a 

1 f+oo fTa -x 2 /i(b 2 +iDt) 

tfj(x,t) = -= / ij(u,t)e iux du 



2vr 7_oo V 2tt y/WTiDi 

It is well known that in this case \ip(x, t)\ 2 has a widening, Gaussian shape all along 
its evolution. We neglect to display pictures of these well known evolutions. 

4.2 Cauchy 

The Cauchy process is one of the most studies non Gaussian, Levy processes [6] , first 
because it is stable, and then because the calculations are relatively accessible. For 
example, if the initial law is a Cauchy with x(u,t) = e~ c *'"', form ( 15 Op and ( T37|) 
we immediately have for the transition chf 

<p(u,t) = e^ +c *)H 

namely the process law remains a Cauchy £b+ct at every t with a typical broadening 
for t — ¥ +oo 

1 b + ct 

P(X,t) = — — rr; r. 18) 

^ ' 7T (b + Ct) 2 + X 2 ^ ' 

Of course this behavior (which is in common with the Gaussian Wiener process) 
comes out from the fact that the Cauchy laws are stable, and we neglect to display 
the corresponding figure. Even when the initial pdf is a Tf,(3) with ipo(u) = (1 + 
&|w|)e~ b ! u ! calculations are easy: now the transition law is again (£ ct , and the one-time 
process law £ ct * T&(3) will have as chf 

t) = X (u, t)<p Q (u) = (1 + b\u\)e~^ u \ 

while the pdf is recovered by chf inversion: 

p(x,t) = — / <b(u,t)e du = — k . (19) 

1 ^ ^J-oo n[(b + ct) 2 + x 2 } 2 K } 

It would be easy to check that this is again a normalized, uni-modal, bell-shaped, 
broadening pdf (see Fig. [3]), with neither an expectation nor a finite variance for 
t > 0. We would find in particular that the process law is the mixture 



of the laws l %$ 1 J 2 {a, 0) of the square root of second kind Beta rw's (see Appendix |E] 

' — 1/2 

for details). Remark that in particular ( B b ' +ct (l/2, 3/2) = T 6+Ct (3). For this example 
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Figure 3: The pdf ffT9l) for a Cauchy process with a Student Tb(3) initial distribution. 



|^r(x,t)| 2 




Figure 4: The square modulus of the Cauchy-Schrddinger wf (12TT) for a Student 
Tft(3) initial distribution. 
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we can also show by direct calculation that the pdfs (1181) and (1191) are both solutions 
of the pseudo-differential Cauchy equation (15TT) . 

The Cauchy-Schrodinger evolutions, on the other hand, show a more interesting 
structure. The simplest case is found when we take as \ipo\ 2 the Student T&(3) 
case P0|) : from fl52|) indeed we have 

^( u ,t) = -f{u,t)Mu) = y/be~ (b+ict)lul 



and hence 



il){x,t) 



ip(u,t)e iux du = \ — 



b + ict 



(21) 



7r (b + ict) 2 + x 2 

This wf (see Figure H]) is correctly normalized in L 2 but shows a new feature: bi- 
modality. In fact \i/j\ 2 has now two well defined maxima smoothly drifting away from 
the center as t — > +oo. It is also possible to show - as an example - that our wf is 
a solution of the Cauchy-Schrodinger equation (j53"|) . For the right-hand side of this 
equation we indeed have from the principal value integral 

ft,, \ , / m c , [2b (b + ict) 2 - x 2 

- [i/j(x + y,t) -^(x,t)\—dy = c\ — 2 

Jy^O V 71 [(b + ICt) 2 + X 2 } 

which is easily seen to coincide with id t ip(x,t). As a consequence the wf ( 12"T|) cor- 
rectly satisfies the pseudo-differential Cauchy-Schrodinger equation flo^]) . A similar 
result is found in the case of a Cauchy Cb initial wf (158"]) : from the propagator (152]) 
we have 

. /oft 

— ir (&M)e- iciH 



^(w,t) 

and hence by inverting the FT: 

il>(x,t) = 



2,7V 

1 



where we defined 



ixybn 



A(z) 



71 



+oo 



.4 



7T 



K {b\u\)e- ictlul e iux du 



x + ct 



A 



x 



ct 



(22) 



| — % arcsinh z 



and we used the following two results 



+oo 



cos(xz)Kn(z) dz 



71 



1 



+oo 



sm(x 'z)Kq(z) dz 



arcsinh x 



x/ 



2 y/l + X 2 

The wf ( |22l is normalized in L 2 and shows (see Figure [5j a behavior similar to that 
of (J21]): its pdf \ip\ 2 starts as a Cauchy distribution and then widens with two 
well defined maxima drifting away from the center. Here too, hence, we have bi- 
modality: remark the difference with the Cauchy process pdfs <£b+ct and C ct * 1^(3) 
which instead broaden by remaining strictly unimodal. 
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0.35 




Figure 5: The square modulus of the Cauchy-Schrodinger wf f|22|) for a Cauchy (£ b 
initial distribution. 



4.3 Laplace 

This bi-modality of the wave packets, or at least its breaking in two symmetric 
structures drifting away from the center can also be found in other examples. Take 
first the Variance-Gamma process of Appendix ID. 31 At variance with the Cauchy 
process, this is an example of a non stable, sd process and hence has a certain 
interest as a non typical case. We will refer to the Appendix IC.4I for a discussion 
of possible initial states. At present we will limit our discussion to initial states of 
the same Variance-Gamma family of the background noise, and we will also always 
choose coincident scale parameters a = b for the background noise and the initial 
states. 

For a Variance-Gamma process with transition law f )54|) and initial pdf fj4Tl) we 
immediately have 



<f>(u,t) = x(u,t)(p (u) 



1 



1 + b 2 u 2 



U+ldt 



and hence the process law simply is 2J(3&(z/ + ut) with pdf 



p{x,t) 



A' 



u+ujt— - 



\X\ 



(23) 



namely always a Variance-Gamma but with a growing parameter v + ut. On the 
one hand this explains why it would be delusory to think of simplifying the example 




Figure 7: The pdf (1231) for a Variance-Gamma process with Variance-Gamma 
53(J5{>(2) initial distribution. 
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-6 
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-2 



Figure 8: The square modulus of the Variance-Gamma-Schrodinger wf (124"]) with a 
Laplace 23(3 6 (1) = £b initial w;/. 



by starting, for instance, with a Laplace = QJ(25&(1) initial law: in fact at every 
time t > the process law would in any case no longer be a Laplace law, but a 
more general Variance-Gamma with u + urt j£ 1. On the other hand this apparently 
explains why at every t the pdf will appear as a broadening, uni-modal distribution 
as shown in the Figures E] and [7] respectively for v — 1 and v — 2. 
For a L-5" evolution, on the other hand, we have from ( I55p and 



i){u,t) 



so that the inverse FT will be 



Y{2u) 



1 



v+iujt 



i r+oo 

Mx,t) = — / ^(u,t)e lux du 

27T 



r(2i/) 



T (2i/ - |) 2 u+iuJt T(u + iut)V2n 
1 / |x| N 

6 It 



a- 



(24) 



Numerical calculations and plotting then show that the wf ( |24"1) always is normalized, 
and that \tp\ 2 has two maxima symmetrically drifting away from the center (see 
Figure |8]). The behavior in x = is rapidly oscillating, but with infinitesimal 



N CUFARO PETRONI: Levy-Schrodinger wave packets 



22 




Figure 9: The square modulus of the Variance-Gamma-Schrodinger wf (|24|) with a 
Variance-Gamma QJ(3b(2) initial wf. 

amplitude as we approach x = 0: in fact the singular behavior of the Bessel function 
is here competing with an infinitesimal \x\ u factor. The distribution shows also a 
slowly decreasing, flat plateau (with micro-oscillations) in the central region, while 
the diverging maxima can be rather dull as in the Figure [9j 

4.4 Poisson 

The following examples come from two ac, but not sd background noises: the com- 
pound Wiener-Poisson processes introduced in the Appendix ID.4I First take the 
process with the transition law $1(2 Dt) * ^3(wt,9T a ) in fl56|) : with a normal initial 
law fl35l) the marginal law of the process becomes yi(2Dt + b 2 ) * *P (out, 9T a ) namely 



OO 



I, .f\k p -x 2 /2{ka?+2Dt+b 2 ) 

p(x, t) = V = (25) 

t^ kl ^2n(ka 2 + 2Dt + b 2 ) 

which apparently is a Poisson mixture of centered, normal pdfs of different vari- 
ances, and hence has the usual bell-like, uni-modal, diffusing shape that we will not 
bother to show. For the other transition law < Xt(2Dt) * ^3 (ut, D a ) in (158]) with the 
same normal initial distribution the mare; inal law instead is m(2Dt + b 2 ) (ut, D a ) 
namely 

~ (ujt)k 1 JL /U\ p -[x-(k-2j)a] 2 /2(2Dt+b 2 ) 
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Figure 10: The pdf fl26|) of a Normal-Poisson process 0t(2.Di) * (ut, D a ) with a 
Gaussian initial law. 



In other words we always have generalized Poisson mixtures, but of non centered 
normal pdf's. Even in this case, however, the shape of the overall pdf will be that 
of a bell-like, uni-modal, diffusing curve (see Figure ITU]) . 

For the L-S equation on the other hand consider first the propagator 9T (2iDt) * 
*P (iut, 9T a ) in fl57j) applied to an initial Gaussian wf fl36l) ; we then have 



ip(u, t) = e 



iuitl e 



^ 2 / 2 -l) 4/262 2+ . m)u2 



TT 



and, by inverting the FT and taking into account the properties of the Gaussian 
integrals, the wf will be 



ij(x,t) = e^J2 



{iuty 



-x 2 /2(ka 2 +2b 2 +2iDt) 



(27) 



, n y/2ir(ka 2 + 2b 2 + 2Wt) 

namely a time-dependent, complex, Poisson superposition of centered Gaussian 
io/'s. The same is true for the second example with propagator 9T (2iDt)*y$ (iut, D a ) 
in ( 1591) with an initial Gaussian wf ( 1361) : the wf FT in fact now is 



i/26 2 _ {b 2 +iDt)u 2 



71 



so that the wf itself will be 



fc=0 



i=o 



J/ ^4^(62 + i£>t) 



(2? 
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Figure 11: The square modulus of the Normal-Poisson Schrodinger wf (128]) with a 
Gaussian initial wf. 

In conclusion, while the plots of p(x, t) in (125]) and ( 126]) simply display the too 
familiar story of a diffusing bell-shaped curve, and the same would be true for 
\ip{x,t)\ 2 in ff271) . for \ip(x,t)\ 2 in (1281) we instead have again a separation of the 
wave packet in two symmetrical sub-packets drifting away from the center (see 
Figure [11]) . 



4.5 Relativist ic qm 

In a way similar to that of the Variance-Gamma, for a Relativistic qm Levy process 
with transition law (160]) and initial distribution (146]) . but with a = 6, we immediately 
have 

4>(u,t) = X (u,t)Mu) = e^< 1 - VT ^ p ) (29) 

p(x,t) = ^ ^ V y 30 

vra v/(z/ + wt) 2 + x 2 /a 2 

and hence the process law simply is %K(v + cut), namely it will stay always in the 
same Relativistic qm family but with a time dependent parameter. The pdf p(x, t) 
is shown in the Figure [12] and has the usual bell-like, uni-modal, diffusing form. 
For the corresponding L-S evolution on the other hand we have from (T4"Tj) and ([6T]) 
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Figure 12: The pdf (|30|) of a Levy process with a Relativistic qm background noise 
and an initial law of the same family. 

that the normalized wf's are 



= iMUu) = y^^y e (^)(i-viw) (31) 

(z/ + iwOe*^ Kl (v'O' + wt) 2 + x 2 /a 2 ) , s 

ih(x t) = , — (32) 

^^(2^) v/(z/ + 2u;t) 2 + x 2 /a 2 

We show in the Figure [T3l how this \i[)(x, y)\ 2 behaves, and in particular, at variance 
with the previous Levy pdf (!30|) . we find here again that the the wf shows two 
symmetric maxima drifting away from the center of the distribution: the bi-modality 
that we have already pointed out in all our other L-S examples. 



5 Conclusions 

We presented in the previous sections several examples of free wave packets that are 
solutions of the L-S equation without potentials (TlTI) . We started by generalizing the 
relation between Brownian motion and Schrodinger equation, and by associating the 
kinetic energy of a physical system to the generator of a symmetric Levy process, 
namely to a pseudo-differential operator whose symbol is the Ich rj of an id law. 
This amounts to suppose, then, that the L-S equation is based on an underlying 
Levy process that can have both Gaussian (continuous) and non Gaussian (jumping) 
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Figure 13: The square modulus of the Relativistic qm wf (132]) with an initial wf of 
the same family. 



components. The use of all the id, even non stable, processes on the other hand 
is important and physically meaningful because there are significant cases that are 
in the domain of our L-S picture, without being in that of the stable (fractional) 
Schrodinger equation. In particular, as discussed in [TJ [6], the simplest form of a 
relativistic, free Schrodinger equation can be associated with a particular type of 
sd, non stable process acting as background noise. Moreover in many instances of 
the Levy-Schrodinger equation the new energy-momentum relations can be seen 
as corrections to the classical relations for small values of certain parameters pQ. 
It must also be remembered that - at variance with the stable, fractional case - 
our model is not tied to the use of processes with infinite variance: the variances 
can be chosen to be finite even in a purely non Gaussian model - as in the case of 
the relativistic, free Schrodinger equation - and can then be used as a legitimate 
measure of the dispersion. Finally let us recall that a typical non stable, Student 
Levy noise seems to be suitable for applications in the models of halo formation in 
intense beam of charged particles in accelerators [H CEO 122] . 

It was then important to explore the general behavior of the diffusing L-S u»/'s: 
we systematically approached this problem by defining in Section [3] a procedure 
allowing us to explore several combinations of initial wf's (Appendix [C]) and back- 
ground Levy noises (Appendix [D]) , and by comparing Levy processes and free L-S 
wave packets. We have then remarked that virtually in all our examples of Section H] 
we witnessed a similar qualitative behavior: first of all the L-S wave packets diffuse, 
in the sense that they broaden in a very regular way. As it is known the variance 
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of a Levy process - when it exists - grows linearly with the time, exactly as in the 
usual diffusions. Of course stable, non Gaussian noises are excluded, since for them 
there is no variance, and we have instead an anomalous sub- and super-diffusive 
behavior. The corresponding L-S wave packets show a similar qualitative behavior 
also if it is not always easy to calculate their variances. 

A second, more surprising feature however is represented by the bi-modality of 
the L-S itf/'s. In fact we found that in virtually all our examples the wave packet 
splits in two sub-packets symmetrically and smoothly drifting away from the center: 
a behavior which is present neither in the free Levy processes, nor in the (Gaussian) 
free Schrodinger wf's. It is interesting to remark, then, that the unique instance 
with a similar bi-modal behavior has been found earlier [TD] deals with confined 
Levy flights. In our opinion the bi-modality found in our examples could then 
be connected to the combined effect of Nelson dynamics, and Levy jumps in the 
background noise, and it would be interesting to explore if this behavior shows up 
again in form of rings and shells respectively for the two- and three-dimensional 
L-S equation. This bi-modality, on the other hand, is in sheer contrast with the 
uni-modality of both the Levy processes and the (Gaussian) Schrodinger wf's. 

It would be important now to explicitly give in full detail the formal association 
between L-S wf's and the underlying Levy processes, namely a true generalized 
stochastic mechanics. In particular we would show that to every wf solution of the 
L-S equation we can associate a well defined Levy process: the techniques of the 
stochastic calculus applied to Levy processes are today in full development [TTJ [T2J 
123] , and at our knowledge there is no apparent, fundamental impediment along this 
road. Finally it would be relevant to explore this Levy-Nelson stochastic mechanics 
by adding suitable potentials to our L-S equation, and by studying the corresponding 
possible stationary and coherent states: all that too will be the subject of future 
papers. 

A Types of laws 

As stated in the Section [1] we deal in this paper with centered laws of ry's X. Even 
when the expectation does not exist we can always speak of centering around the 
median. On the other hand to eliminate the centering it will be enough to take 
X + b with 6 e R instead of X, then to substitute x — b to x in the /, and to add 
a factor e lbu to the chf ip. For our purposes it will also be expedient to introduce a 
dimensional scale parameter a > to take into account the physical dimensions of 
our rv's: to fix the ideas in this paper a will be supposed to be a length. Take first a 
tv X with law pdf f and chf <p, and suppose that X is a dimensionless quantity; 
then the variables argument of / and (p, will be dimensionless. On the other hand 
X a = aX will be a length and will follow a law 3a with 
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Here x and u are now dimensional variables (x is a length, while u is the reciprocal 
of a length), so that x/a and au will be dimensionless. Remark that within this 
notation we numerically have $ = $1, s° that for instance fi(x) = f(x). This 
could be slightly misleading since the argument of f\ is a length, while that of / 
is supposed to be dimensionless. To avoid any possible misunderstanding we will 
then reserve the symbols / and (p for the dimensionless laws, while fi and (p± 
will be associated to the dimensional ones. For example if X follows the standard, 
dimensionless normal law 9T with 

f{x) = 6 —=- , <p( u ) = e^ 2 



'2-n 

the dimensional rv's X a = aX will follow the laws OTq, = 9T(a 2 ) with 

—x 2 /2a? 

f a {x) = e —j= , <p a (u) = e- a2u2 l\ 



ay2n 

Then / and fx will be coincident, but the dimensional meaning of their respective 
variables will be different. Remark finally that in general we will choose dimension- 
less laws that are not necessarily standard laws: of course (when the variances exist) 
we will have V[X ] = a 2 V[X], but V[X] is not always supposed to be equal to 1. 

We could now think to $ a as the parametric family of the rescaled rw's aX: 
these parametric families spanned just by one scale parameter a are here entire 
types of kmsB in fact, since here we only deal with centered laws (see Section [TJ, no 
centering parameter b is required, and our types are spanned by means of the scale 
parameter a only. In this paper we will also consider other parametric families of 
laws with some dimensionless parameter A, which will not in general be coincident 
with the scale parameter a. We could then have two-parameters families 3a(A), 
and in general we are interested in finding which sets are closed under convolution 
(namely under addition of the corresponding independent rv's). When a type of 
laws is closed under convolution (as in the normal case of the previous example) its 
laws are said to be stable: the convolution would produce another law of the same 
type, namely a law with only a different scale parameter (in our notation: same A, 
but different a). If instead the convolution produces a law of the same family, but 
not of the same type (different A), then the family is closed under convolution, but 
its laws are not stable: this is the case, among others, of the Variance-Gamma laws 
QJ(3 a (A). Finally, when the result of a convolution is a law not belonging at all to 
the family, then 5a(A) is not even closed under convolution, as for the Student X a (A) 
family. 



4 A type of laws (see [16] Section 14) is a family of laws that only differ among themselves by a 
centering and a rescaling: in other words, if tp(u) is the chf of a law, all the laws of the same type 
have c/i/'s e zhu (p(au) with a centering parameter b £ M, and a scaling parameter a > (we exclude 
here the sign inversions). In terms of ry's this means that the laws of X and aX + b (for a > 0, 
and b £ R) always are of the same type, and on the other hand that X and Y belong to the same 
type if and only if it is possible to find a > 0, and b £ R such that Y and aX + b have the same 
law, namely Y = aX + b. 
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B Symmetric and ac, compound Poisson laws 

Among the id, non sd laws the Poisson case stands as the most important example, 
but the simple Poisson law is neither symmetric, nor ac. We will then generalize 
it in order to avoid these shortcomings. A Poisson law ^J(A) is a non symmetric, 
non sd, non ac, id law without Gaussian component (/3 = 0). The probability is 
concentrated on the integer numbers with the usual Poisson distribution so that 
formally 

00 \ k 

f(x) = h(x) , tp(u) = eH^-V , £(x) = A 5 1 (x) . 

fc=0 

Both expectation and variance have value A. Since ^J(A) is neither centered, nor 
symmetric the generator of the corresponding Levy process will not be self-adjoint. 
It is well known, moreover, that the sample paths of the corresponding simple Pois- 
son process are ascending staircase trajectories, with randomly located steps of unit 
height, A representing the average number of jumps per unit time interval. As a con- 
sequence these processes are not ac. To move ahead we must then first symmetrize 
the Poisson law, and then make it ac. 

Take a symmetric (we do not require it to be ac or id) law 9) with chf = 
and build the corresponding compound Poisson law ^3(A, Sji) with chf 

v (u) = e - x V ^ #*(u) = e mu) ~ 1] , rj(u) = A[#(u) - 1] 

k=0 

thus generalizing the simple Poisson case where ~&(u) = e m . When S) is also ac with 
pdf h(x) the law of <p(A, Sj) is 

{k times 
h*...*h, k = l,2... (33) 
S , k = 

but we can immediately see that this is still not ac even if has a density: in fact 
for k = we always have a degenerate law 5 . The compound Poisson law ^P(A, S)) 
has neither a drift (a = because of the required symmetry) nor a Gaussian part 
(j3 — 0), and its Levy pdf (that we will suppose for simplicity to show no singularities 
at x = 0) is £(x) = Xh(x): namely we have C = (0, 0, Xh). The laws of the increments 
of the corresponding compound Poisson process ^(ujt,Sj) with u> = A/r are then 
the time dependent mixtures 



00 



p( x ,t) = c- t J2 i ^fh* k (x) (34) 

k=0 



while its self-adjoint generator (no singularities are present at x — 0) is 

/+00 
[v(x + y) — v(x)\h(y) dy . 
■00 
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Its sample trajectories are now up and down staircase functions, with steps at Pois- 
son random times, and random jump heights distributed according to the symmetric 
law .fj. Since however for k = the law is degenerate in x = 0, these sample tra- 
jectories stick at x = for a finite time (with probability 1), and the marginal 
distribution of the process is not ac. In other Levy processes instead (as the Wiener 
process for example) the trajectory starts at x — 0, but its random path immediately 
leaves this position. 



To give a first example of these symmetric (but not ac) compound Poisson laws 
take S) = Vl a so that h* k ~ 0T(/ca 2 ) for k — 0, 1, . . . ; we then have for <p(A, 9t a ) 



The transition pdfs of the corresponding compound Poisson process are then the 
time dependent mixtures of yi(ka 2 ) laws 



As another example suppose instead that $) = D a is a Bernoulli symmetric law, 
doubly degenerate around the positions ±a, namely 





and the generator takes the form 




h(x) 




+ 5- 



-i 




~&(u) = cos au, 



and remark that now 




As a consequence we will have for *P(A, D a ): 




r)(u) 
£(x) 



A(cos au — 1) 




We then easily have for the transition law of the process 
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while the generator is 

[Av](x) = ^ [v(x + 1) - 2v(x) + v(x - 1)] 

We will then further generalize our compound Poisson distributions in order to 
get ac laws and processes. Take a compound Poisson law ^(X,^), and another 
independent, symmetric, ac, id law $) with pdf h (x), chf $o( u ) — c^ ^ and Levy 
triplet £ = (0,(3 ,£ ). Consider then the law i) * ^P(A,i}) obtained by addition 
(convolution) so that 

<f(u) = M^)e mu) - 1} , V(u) = Co(«) + A(0(u) - 1) 

while the pd/ is 

{fc times 
h *h7^~7h, k = 1,2... 
/i , = 

This is now a mixture of ac laws. The law Sjq * ^3(A, Sj) will also be symmetric if 
both h and /io are symmetric, and it will have a Gaussian component if /3o 7^ 0. As 
a consequence we will have a = from the symmetry, £(x) = Xh(x) + £q(x), and 
finally £ = (0 , O , Xh + £ ). The laws of the increments of the corresponding Levy 
process will then be ip(t) = e xt ^~ 1 ^ T , namely 

770, t) = - Co(«) + - 1] 

r 

so that the process will be the superposition of two independent processes: an S) - 
Levy process plus a ^P(ut, Sj) compound Poisson process. Its trajectories will then 
be the paths of the f)o~Levy process, interspersed with Poisson random jumps with 
size law ij. If then h (x, t) is the pdf of $q T (u), the t-increment pdfs of our process 
will be 

k\ 



00 



p( x ,t) = e-" t Y,^n- [ho(t)*h* k ] (x) 



k=0 

and the self-adjoint process generator 

[Av](x) = %.%v{x)+ [ [v(x + y)-v(x)][Xh(y) + e (y)}dy. 



2 



Possible examples of these Sjq are both the Gaussian and the non Gaussian stable 
laws (in particular the Cauchy process), and several self-decomposable laws as the 
Student or the Variance-Gamma. The relevant particular case of a Gaussian Sjo is 
discussed in the Section 12.41 
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C Initial states 

We define here a list of possible initial pdf's and io/'s. To simplify our calculations we 
will choose the initial pdfs to be centered and symmetric, and whenever convenient 
we will take pairs fo, ipo satisfying the relation f = \ipo\ 2 . Remark that, while fo is a 
normalized (in L 1 ) pdf and ipo is a (non normalized, and possibly non normalizable) 
chf ip and ip must be both normalized (in L 2 ) wf's so that we must always pay 
attention to the constants which are in front of them. Here moreover - to put in 
evidence the meaning of the involved quantities - our laws and time coordinates will 
be dimensional: the space a, b and time r scaling parameters will be explicitly taken 
into account. 



C.l Normal 91 & 

Initial laws and wf's with fo = \ipo\ 2 are in this case 

p -x 2 /2b 2 

fo(x) = -j=, = e- fe2 " 2/2 (35) 

-x 2 /4b 2 /oL2 

M*) = ^TTI ' Mu) = ti— e- fe2 " 2 (36) 

V27r6 2 V 71 

Remark that, while ip is just the square root of f , ip is the FT of ip an d its 
relation to <p is given by the equation (JTS]) . The two wfs, moreover, are both 
normalized in L 2 . 



C.2 Cauchy U = % b (l) 
Initial laws and wfs in this case are 



fo(x) = -J- -^— 2 , Mu) = (37) 

07T Z + X z 



$o(x) = -fi=J IT"! — a ' = ~~ K o(b\u\) (31 

Vbn V b z + x z 7i 



where K is the modified Bessel function of order 0; it is easy to show indeed that 
(see for example [20] 9.6.21) 



b 2 \/2b 
^ Q (u) = ^=-^= / \l-r 2 ;e- t » ai dx = —K {b\u\). 



The normalization H^oll 2 = lj and the relation (po = ipo * ipo, are then 

K 2 (b\u\)du = ^ 



K (b\u-v\)K (b\v\)dv = —e 



TT 2 
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The first can be reduced to 

r+oo 2 
J o K 2 (u)du = - 

which can be verified by direct calculation. On the other hand the convolution, that 
can be reduced to the dimensionless relation 



/ K (\u-v\)K (\v\)dv = ^e-\ u \, 



does not seem to be an otherwise known result. 

C.3 3-Student T 6 (3) 
Initial laws and iw/'s in this case are 



fo(x) = I (^) 2 , Mu) = e- b ^(l + b\u\) (39) 

Ms) = \/r- UTT » Mu) = Vbe'W (40) 



bn b 2 + x 2 ' 

It is very easy to show that is the right FT of ipo 

^b(u) = J-— —— e— dx = Vbe^ 
V bn y/2-K J-oo b 2 + x 2 

while here again an elementary calculation shows also that (po = tpo * i>o- 
C.4 Variance— Gamma QJ(5b(z/) 

In the general Variance-Gamma case, to make calculations possible, we will not 
always choose pairs of initial pdf's and io/'s satisfying ifj = vTo- A possible example 
then is 



\x\ \ z „ / UC 



/oW = n ,.,, ^ 1 V) **-i(V)> ^)=(r^o) (4i) 



2»r(v)V2^b\b J "-*\b J' ^ uv 7 \l + b 2 u< 

^-/v^iy(™)" (43) 

where the functions are chosen in order to have an evolution easy to calculate. The 
wf's ipo and i/jq, in any case, are both normalized in L 2 (as can be seen by direct 



N CUFARO PETRONI: Levy-Schrodinger wave packets 



34 



calculation) and are apparently in the FT relation. As a consequence here the Levy 
and the L-S evolutions will possibly start with different pdf's. In fact the usual 
relation f = \ipo\ 2 could be easily restored just in the particular case of v — 1, 
namely for an initial Laplace law £{, = 03(25 b(l): 

e -M/& i 

AM = — , Mu) = Y+Hhf (44) 

e -M/2& „ 2 

*> W = W * (M) = VxTT4W («) 

Here it is elementary to check indeed that ip$ = a/Jo, that is the FT of ?/>o> an d 
finally that — * V'o- This particular case, however, is not really easier than 
the general case of the Variance-Gamma process. In fact, as we will see soon, the 
parameter affected by the time evolution is exactly u, so that it is of no help to start 
with v — 1 if it immediately becomes v ^ 1. 

C.5 Relativistic qm £H&(^) 

Again to make calculations easy we will choose as initial chf and wf FT respectively 



ve v K x ^v 2 + x 2 /b 2 



fo(x) = - y = ~A = e^ 1 -^/^) (46) 

Ux) = ^ptA ' ■ o(,i) = V ™ e ^ (47) 

which are in a relation similar to that of ( [41]) - ( 143]) . It is easy to recognize that the 
wf's are correctly normalized in L? . 

D Transition laws and propagators 

We will list here a few examples of background Levy noises by paying attention 
to pick up processes with a known transition pdf associated to the evolution equa- 
tion (|T5|) and a known propagator associated to the free L-S equation ffTTj) . 

D.l Normal 9T(2L>£) 



Here the background noise is a Wiener process: take a DT a law with Levy triplet 
£ = (0,a,0) 



-x 2 /2a 2 
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The transition law of the corresponding Levy process is then $1(2 Dt) with D ■ 
a 2 /2r, namely 

-x 2 /ADt 

q(x,t) = -j=, X (u,t) = e~ Dtu2 (4f 



and the pdf evolution equation ([13]) is the usual Fokker-Planck equation 

d t p(x, t) = D dlp(x, t) 

The corresponding L-S propagator 0X(2iZ)t) is again formally normal albeit with an 
imaginary variance: 

-x 2 /4iDt 

g{x,t) = -j=, 7(M) = e- !Dt " 2 (49) 

and hence the L-S equation (|T7j) is the usual free Schrodinger equation 

id t ip(x,t) = —Ddlip(x,t). 



D.2 Cauchy £ 



ct 



From the Cauchy law <£ a , a typical stable, non Gaussian law with Levy triplet 
C = (0, 0, a/ ttx 2 ) and with 

M = - -^—^ ¥>(«) = e- a|w| 

an cr + ar 

we get the transition law € ct of the Cauchy process with c = ajr: 

<l^t) = —-^— 2l X (u,t) = e-^ (50) 

TTCl C z t z + X z 

and the corresponding process equation ( TT3l 

d t p{x,t) = \p(x + y,t) - p(x,t)]-^—dy. (51) 
•Vo Try 

On the other hand the L-S propagator £ ict is 

fl(M) = - 2 , 7(«,t) = e- ldW (52) 

Z7T trr z — x z 



and the L-5" equation (JT 



c 



id t i/j(x,t) = - / [V>(a: + y,*) - ^(x,t)] — rdy (53) 
'^o Try 



Remark that, at variance with the transition pdf (1501) . the Cauchy-Schrodinger 
propagator (152|) has two simple poles in x = ±ct drifting away from the center 
x = with velocity c. 
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D.3 Variance— Gamma %3(5 a (ujt) 

Take a sd, non stable Variance-Gamma law QJ(J5 a (A) with symmetric Levy triplet 
£ = (0, 0, Xe-^ a /\x\) and with 

/(*) = 7^ f-V V± (- ) ■ 



2 A r(A) v / 2^a V a J 5 V a /' Vi + ^w 2 

The transition law will then be QJC5 a (a;t) with uj = A/r: 

?M = (t)"' (t) • * M) ^ (ot)"* < 54 » 

and the corresponding process equation ([TBI) 

/■ e -|tfl/« 
d t p(x, t)=u \p{x + y,t)- p{x, t)] — j— — . 

Jy^o \y\ 

so that the evolution will only affect the parameter A, while a will always be the 
same. Then for the L-S propagator ^XS(& a {iujt) we have 

I i\ iujt— i s\ i\ / 1 \ 

# \ / kc \ , / l 



2 4aJ *r(?cc;t)\/27ra \ a / 2 \ a / y 1 + a u 

(55) 

while the L-S 1 equation (ITT)) becomes 

/■ e -M/ a 
id t ip{x, t) = —u I [ip(x + y,t) - ip(x, £)] — — — rfy 

</^o \y\ 

D.4 Wiener-Poisson 9I(2D£) * (otf, #) 

We will consider here two examples of id, non sd background noise: for notations and 
details see Section EH and Appendix [Bl Take first the law 9to. * (A, 9T a ) discussed 
in the Section 12.41 From its chf we see that, with oj = A/r and D = cr 2 /2r, the 
transition law W(2Dt) * (ut, 9? a ) 

°° f, „-x 2 /2(ka 2 +2Dt) , 2 

g (x,t) = e-^y;M- g x(«,*)=e^ e -"" ^e"^ 2 (56) 

The corresponding Wiener-Poisson process will have sample paths which are Brow- 
nian trajectories interspersed with Gaussian jumps at Poisson times with intensity 
A. The process pdf's then have an elementary form as time dependent Poisson mix- 
tures of time dependent normal laws and the corresponding process equation ( TT3|) 
will become 

/+oo e -y 2 /2a 2 
\p(x +y,t)- p(x, t)] dy. 
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The L-S propagator 9t (2iDt) * (icot, 9T a ) now is 

°° (ii iA k -x 2 /2(ka 2 +2iDt) „„, 

and the L-5 equation ( IT?]) becomes 



id t iP(x,t) = -Dd 2 x iP(x,t) - u / [V(g + y,t)--0(g,t)] / - — - dy. 

As a second example take the law 9^ * (A, D a ) discussed in the Section I2.4[ from 
its Ich rj(u, t) = ut(cos au — 1)— Dtu 2 we see that the law of the corresponding Levy 
process is < fl(2Dt) * ^(out, D a ) and hence 

X(u,t) = e ^(cosa U -l)-OtM 2 

while the process equation (IT31 is 



o * nA 2 / #\ . p(x + a,t)-2p(x,t)+p(x-a) 
d t p{x, t) = Dd x p(x, t) + u . 

This process will have sample paths which are again Brownian trajectories inter- 
spersed with jumps ±0 at Poisson times with intensity A. The L-S propagator 
m(2iDt) * %$(iut, O a ) instead is 

^ (iut) k 1 Jt. fk\ e -[x-(k-2j)a] 2 /4iDt 
ry( u t) = gi^ticosau-l^iDtu 2 

and the L-S equation ( 1T7|) is 

id t i>(x,t) = —Ddlip(x,t) —00 



, ; ip(x + a, t) — 2ip(x, t) + ip(x — a) 



D.5 Relativistic qm 9K a (ujt) 

We immediately see from the chf of 9^ a (A) that the corresponding Levy process 
9\ a (u)t) will have as transition law 



ute^Ki ^u 2 t 2 + x 2 /a 2 



t) = ) 2+2 : =. = \ X(u, t) = e^-^ 1 ^) (60) 

ixa^JuH* + x 1 jar 
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with u = X/t as usual. We can also explicitly write the process equation ( Tl3~ 
d t p{x,t)=u [ [p(x + y,t) -p(x,t)] Kl f a dy. 

Jy^o n\y\ 

On the other hand the L-S propagator 9\ a (iujt) will be given by 



9(x,t) = / 2+2 A l{u,t) = eKi-v^^) (61) 

ixay —u z t z + x z /a z 

with singularities in x = ±aut, and corresponds to the L-S equation 

id t ^(x,t) = -u [ &(x + y,t)-iP(x,t)} Kl{l f( a) dy. 

Jy^o n\y\ 



We remember here, as remarked in the Section 12. 5[ that this essentially is the 
integro-differential form of the well known relativistic, free Schrodinger equation 

ihd t ijj(x,t) = i/m 2 c 4 — c 2 h 2 d 2 ip{x,t) 

that we recover by taking u> = X/t = mc 2 /h, a = h/mc, and by reabsorbing an 
irrelevant constant term mc 2 in a phase factor of the wf [1] . 

E Second kind Beta laws 

If Z is a rv with a (dimensionless) Beta law 05 (a, (3) (a, j3 > 0) namely with pdf 

fz(z) = ^ ( , < z < 1 

then Y = Z/(l — Z) is distributed according to a second kind Beta law 55 (a, /?) with 
(see for example [21]) 

/y(l/) = 5(ab)(ll^' 

We could also introduce a scale parameter a to get the types 23 a (a, /3) and 25 a (a, /3), 
but to simplify the notation we will first consider only the dimensionless laws. Take 
now a third rv X = e\/Y where \/Y is the positive square root of Y, while e is 
another independent rv taking the two values ±1 with the same probability 1/2. 
We find then that its pdf is 



fx(x) 



(X 2 ) a - L 2 



B(a,P) (l + x 2 ) 



a+/3 
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We will use for these laws the symbol !B 1 / 2 (a, (5) because X is the square root of a 
second kind Beta rv. In particular we recover the family of the Student laws as 

® 1/2 Q^)=*(A), A >°> 

while Q3 1 / 2 (3/2 , 1/2) is the law introduced in the Section l4"72l to describe the evolution 
of an initial Student law T(3) by a Cauchy transition pdf. For this law we have 

= - 77-^, = (l-M)e-W, 

TV (1 + X z ) z 

and (as for the usual Cauchy laws) we find that it has neither an expectation, nor a fi- 
nite variance. The decomposition ( 12"U|) could now be written also as a relation within 
the (dimensional) family 2$y 2 (a,/3) by remembering that £ ct = s B 1 c [ 2 (l/2 , 1/2) and 

— 1/2 

% = <8l'*(l/2, 3/2). In fact, for given arbitrary scale parameters a and b, and with 



2 a + b 

W 2 a + 6 2 \ a + b J 7 
we easily see that (1201) is a special case (for a = ct) of 

»^ 6- 1) • ^ (I- 1) = (|. I) + « (|. I) 

namely 

(1 + 6| M |)e- (a+fe)|n| = P [1 - (a + &)M]e-( a+b )l"l + Q [1 + (a + 6)| M |] e - (a+fe)|n| 
(a + 6) 2 (a + 26) + ax 2 2P (a + b)x 2 2Q (a + b) 3 
7r[(a + 6) 2 + a; 2 ] 2 ~ V [(a + b) 2 + x 2 } 2 + ~V [{a + b) 2 + x 2 } 2 
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